Human cytomegalovirus extensively re-organizes the human genome, diminishing TEAD1 transcription factor activity.

Human cytomegalovirus (HCMV) infects up to 80% of the world's population. Here, we show that HCMV infection leads to widespread changes in human chromatin accessibility and chromatin looping, with hundreds of thousands of genomic regions affected 48 hours after infection. Integrative analyses reveal HCMV-induced perturbation of Hippo signaling through drastic reduction of TEAD1 transcription factor activity. We confirm extensive concordant loss of TEAD1 binding, active H3K27ac histone marks, and chromatin looping interactions upon infection. Our data position TEAD1 at the top of a hierarchy involving multiple altered important developmental pathways. HCMV infection reduces TEAD1 activity through four distinct mechanisms: closing of TEAD1-bound chromatin, reduction of YAP1 and phosphorylated YAP1 levels, reduction of TEAD1 transcript and protein levels, and alteration of TEAD1 exon-6 usage. Altered TEAD1-based mechanisms are highly enriched at genetic risk loci associated with eye and ear development, providing mechanistic insight into HCMV's established roles in these processes.

Shared and distinct interactions of type 1 and type 2 Epstein-Barr Nuclear Antigen 2 with the human genome.

BACKGROUND: There are two major genetic types of Epstein-Barr Virus (EBV): type 1 (EBV-1) and type 2 (EBV-2). EBV functions by manipulating gene expression in host B cells, using virus-encoded gene regulatory proteins including Epstein-Barr Nuclear Antigen 2 (EBNA2). While type 1 EBNA2 is known to interact with human transcription factors (hTFs) such as RBPJ, EBF1, and SPI1 (PU.1), type 2 EBNA2 shares only ~ 50% amino acid identity with type 1 and thus may have distinct binding partners, human genome binding locations, and functions. RESULTS: In this study, we examined genome-wide EBNA2 binding in EBV-1 and EBV-2 transformed human B cells to identify shared and unique EBNA2 interactions with the human genome, revealing thousands of type-specific EBNA2 ChIP-seq peaks. Computational predictions based on hTF motifs and subsequent ChIP-seq experiments revealed that both type 1 and 2 EBNA2 co-occupy the genome with SPI1 and AP-1 (BATF and JUNB) hTFs. However, type 1 EBNA2 showed preferential co-occupancy with EBF1, and type 2 EBNA2 preferred RBPJ. These differences in hTF co-occupancy revealed possible mechanisms underlying type-specific gene expression of known EBNA2 human target genes: MYC (shared), CXCR7 (type 1 specific), and CD21 (type 2 specific). Both type 1 and 2 EBNA2 binding events were enriched at systemic lupus erythematosus (SLE) and multiple sclerosis (MS) risk loci, while primary biliary cholangitis (PBC) risk loci were specifically enriched for type 2 peaks. CONCLUSIONS: This study reveals extensive type-specific EBNA2 interactions with the human genome, possible differences in EBNA2 interaction partners, and a possible new role for type 2 EBNA2 in autoimmune disorders. Our results highlight the importance of considering EBV type in the control of human gene expression and disease-related investigations.

Association between systemic lupus erythematosus and myasthenia gravis: A population-based National Study.

BACKGROUND: Systemic lupus erythematosus (SLE) and myasthenia gravis (MG) are autoimmune diseases. Previous case reports and case series suggest an association may exist between these diseases, as well as an increased risk of SLE after thymectomy for MG. We undertook this study to determine whether SLE and MG were associated in large cohorts. METHODS: We searched the IBM Watson Health Explorys platform and the Department of Veterans Affairs Million Veteran Program (MVP) database for diagnoses of SLE and MG. In addition, we examined subjects enrolled in the Lupus Family Registry and Repository (LFRR) as well as controls for a diagnosis of MG. RESULTS: Among 59,780,210 individuals captured in Explorys, there were 25,750 with MG and 65,370 with SLE. 370 subjects had both. Those with MG were >10 times more likely to have SLE than those without MG. Those with both diseases were more likely to be women, African American, and at a younger age than MG subjects without SLE. In addition, the MG patients who underwent thymectomy had an increased risk of SLE compared to MG patients who had not undergone thymectomy (OR 3.11, 95% CI: 2.12 to 4.55). Autoimmune diseases such as pernicious anemia and miscellaneous comorbidities such as chronic kidney disease were significantly more common in MG patients who developed SLE. In the MVP, SLE and MG were also significantly associated. Association of SLE and MG in a large SLE cohort with rigorous SLE classification confirmed the association of SLE with MG at a similar level. CONCLUSION: While the number of patients with both MG and SLE is small, SLE and MG are strongly associated together in very large databases and a large SLE cohort.

Lupus enhancer risk variant causes dysregulation of IRF8 through cooperative lncRNA and DNA methylation machinery.

Despite strong evidence that human genetic variants affect the expression of many key transcription factors involved in autoimmune diseases, establishing biological links between non-coding risk variants and the gene targets they regulate remains a considerable challenge. Here, we combine genetic, epigenomic, and CRISPR activation approaches to screen for functional variants that regulate IRF8 expression. We demonstrate that the locus containing rs2280381 is a cell-type-specific enhancer for IRF8 that spatially interacts with the IRF8 promoter. Further, rs2280381 mediates IRF8 expression through enhancer RNA AC092723.1, which recruits TET1 to the IRF8 promoter regulating IRF8 expression by affecting methylation levels. The alleles of rs2280381 modulate PU.1 binding and chromatin state to regulate AC092723.1 and IRF8 expression differentially. Our work illustrates an integrative strategy to define functional genetic variants that regulate the expression of critical genes in autoimmune diseases and decipher the mechanisms underlying the dysregulation of IRF8 expression mediated by lupus risk variants.

A High Prevalence of Anti-EBNA1 Heteroantibodies in Systemic Lupus Erythematosus (SLE) Supports Anti-EBNA1 as an Origin for SLE Autoantibodies.

BACKGROUND: That Epstein-Barr virus (EBV) infection is associated with systemic lupus erythematosus (SLE) is established. The challenge is to explain mechanistic roles EBV has in SLE pathogenesis. Previous studies identify four examples of autoantibody cross-reactions between SLE autoantigens and Epstein-Barr nuclear antigen 1 (EBNA1). For two of these examples, the earliest detected autoantibody specifically cross-reacts with EBNA1; thereby, defined EBNA1 epitopes induce a robust autoantibody response in animals. These results suggest that the autoantibodies initiating the process leading to SLE may emerge from the anti-EBNA1 heteroimmune response. If this hypothesis is true, then anti-EBNA1 responses would be more frequent in EBV-infected SLE patients than in EBV-infected controls. We tested this prediction. METHODS: We evaluated published East Asian data by selecting those with a positive anti-viral capsid antigen (VCA) antibody immunoglobulin G (IgG) test and determining whether anti-EBNA1 was more common among the EBV-infected SLE cases than among matched EBV-infected controls with conditional logistic regression analysis. RESULTS: All the qualifying SLE patients (100%) in this dataset were EBV-infected compared to age- and sex-matched controls (92.2%) [odds ratio (OR) = 28.6, 95% CI 6.4-∞, p = 8.83 × 10-8], confirming the known close association of EBV infection with SLE. Furthermore, virtually all the SLE cases have both anti-VCA IgG and anti-EBNA1 IgG antibodies [124 of 125 (99.2%)], which are more frequently present than in age- and sex-matched EBV-infected controls [232 of 250 (93.2%)] (OR = 9.7, 95% CI 1.5-414, p = 0.0078) for an 89.7% SLE attributable risk from anti-EBNA1, which is in addition to the 100% SLE risk attributable to EBV infection in these data. CONCLUSIONS: The association of EBV infection with SLE is reconfirmed. The prediction that anti-EBNA1 is more frequent in these SLE cases than in EBV-infected controls is true, consistent with the hypothesis that anti-EBNA1 contributes to SLE. This second EBV-dependent risk factor is consistent with a molecular mimicry model for the generation of SLE, starting with EBV infection, progressing to anti-EBNA1 response; then molecular mimicry leads to anti-EBNA1 antibodies cross-reacting with an SLE autoantigen, causing autoantibody epitope spreading, and culminating in clinical SLE. These results support the anti-EBNA1 heteroimmune response being a foundation from which pathogenic SLE autoimmunity emerges.

The US Department of Veterans Affairs Science and Health Initiative to Combat Infectious and Emerging Life-Threatening Diseases (VA SHIELD): A Biorepository Addressing National Health Threats.

Background: The coronavirus disease 2019 (COVID-19) pandemic, caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has demonstrated the need to share data and biospecimens broadly to optimize clinical outcomes for US military Veterans. Methods: In response, the Veterans Health Administration established VA SHIELD (Science and Health Initiative to Combat Infectious and Emerging Life-threatening Diseases), a comprehensive biorepository of specimens and clinical data from affected Veterans to advance research and public health surveillance and to improve diagnostic and therapeutic capabilities. Results: VA SHIELD now comprises 12 sites collecting de-identified biospecimens from US Veterans affected by SARS-CoV-2. In addition, 2 biorepository sites, a data processing center, and a coordinating center have been established under the direction of the Veterans Affairs Office of Research and Development. Phase 1 of VA SHIELD comprises 34 157 samples. Of these, 83.8% had positive tests for SARS-CoV-2, with the remainder serving as contemporaneous controls. The samples include nasopharyngeal swabs (57.9%), plasma (27.9%), and sera (12.5%). The associated clinical and demographic information available permits the evaluation of biological data in the context of patient demographics, clinical experience and management, vaccinations, and comorbidities. Conclusions: VA SHIELD is representative of US national diversity with a significant potential to impact national healthcare. VA SHIELD will support future projects designed to better understand SARS-CoV-2 and other emergent healthcare crises. To the extent possible, VA SHIELD will facilitate the discovery of diagnostics and therapeutics intended to diminish COVID-19 morbidity and mortality and to reduce the impact of new emerging threats to the health of US Veterans and populations worldwide.

Desmoplakin and periplakin genetically and functionally contribute to eosinophilic esophagitis.

Eosinophilic esophagitis (EoE) is a chronic allergic inflammatory disease with a complex underlying genetic etiology. Herein, we conduct whole-exome sequencing of a multigeneration EoE pedigree (discovery set) and 61 additional multiplex families with EoE (replication set). A series of rare, heterozygous, missense variants are identified in the genes encoding the desmosome-associated proteins DSP and PPL in 21% of the multiplex families. Esophageal biopsies from patients with these variants retain dilated intercellular spaces and decrease DSP and PPL expression even during disease remission. These variants affect barrier integrity, cell motility and RhoGTPase activity in esophageal epithelial cells and have increased susceptibility to calpain-14-mediated degradation. An acquired loss of esophageal DSP and PPL is present in non-familial EoE. Taken together, herein, we uncover a pathogenic role for desmosomal dysfunction in EoE, providing a deeper mechanistic understanding of tissue-specific allergic responses.

Validation of low-coverage whole-genome sequencing for mitochondrial DNA variants suggests mitochondrial DNA as a genetic cause of preterm birth.

Preterm birth (PTB), or birth that occurs earlier than 37 weeks of gestational age, is a major contributor to infant mortality and neonatal hospitalization. Mutations in the mitochondrial genome (mtDNA) have been linked to various rare mitochondrial disorders and may be a contributing factor in PTB given that maternal genetic factors have been strongly linked to PTB. However, to date, no study has found a conclusive connection between a particular mtDNA variant and PTB. Given the high mtDNA copy number per cell, an automated pipeline was developed for detecting mtDNA variants using low-coverage whole-genome sequencing (lcWGS) data. The pipeline was first validated against samples of known heteroplasmy, and then applied to 929 samples from a PTB cohort from diverse ethnic backgrounds with an average gestational age of 27.18 weeks (range: 21-30). Our new pipeline successfully identified haplogroups and a large number of mtDNA variants in this large PTB cohort, including 8 samples carrying known pathogenic variants and 47 samples carrying rare mtDNA variants. These results confirm that lcWGS can be utilized to reliably identify mtDNA variants. These mtDNA variants may make a contribution toward preterm birth in a small proportion of live births.

Complement genes contribute sex-biased vulnerability in diverse disorders.

Many common illnesses, for reasons that have not been identified, differentially affect men and women. For instance, the autoimmune diseases systemic lupus erythematosus (SLE) and Sjögren's syndrome affect nine times more women than men1, whereas schizophrenia affects men with greater frequency and severity relative to women2. All three illnesses have their strongest common genetic associations in the major histocompatibility complex (MHC) locus, an association that in SLE and Sjögren's syndrome has long been thought to arise from alleles of the human leukocyte antigen (HLA) genes at that locus3-6. Here we show that variation of the complement component 4 (C4) genes C4A and C4B, which are also at the MHC locus and have been linked to increased risk for schizophrenia7, generates 7-fold variation in risk for SLE and 16-fold variation in risk for Sjögren's syndrome among individuals with common C4 genotypes, with C4A protecting more strongly than C4B in both illnesses. The same alleles that increase risk for schizophrenia greatly reduce risk for SLE and Sjögren's syndrome. In all three illnesses, C4 alleles act more strongly in men than in women: common combinations of C4A and C4B generated 14-fold variation in risk for SLE, 31-fold variation in risk for Sjögren's syndrome, and 1.7-fold variation in schizophrenia risk among men (versus 6-fold, 15-fold and 1.26-fold variation in risk among women, respectively). At a protein level, both C4 and its effector C3 were present at higher levels in cerebrospinal fluid and plasma8,9 in men than in women among adults aged between 20 and 50 years, corresponding to the ages of differential disease vulnerability. Sex differences in complement protein levels may help to explain the more potent effects of C4 alleles in men, women's greater risk of SLE and Sjögren's syndrome and men's greater vulnerability to schizophrenia. These results implicate the complement system as a source of sexual dimorphism in vulnerability to diverse illnesses.

Pleiotropy in the Genetic Predisposition to Rheumatoid Arthritis: A Phenome-Wide Association Study and Inverse Variance-Weighted Meta-Analysis.

OBJECTIVE: This study was undertaken to investigate the hypothesis that a genetic predisposition toward rheumatoid arthritis (RA) increases the risk of 10 cardiometabolic and autoimmune disorders previously associated with RA in epidemiologic studies, and to define new genetic pleiotropy present in RA. METHODS: Two approaches were used to test our hypothesis. First, we constructed a weighted genetic risk score (wGRS) and then examined its association with 10 prespecified disorders. Additionally, a phenome-wide association study (PheWAS) was carried out to identify potential new associations. Second, inverse variance-weighted regression (IVWR) meta-analysis was used to characterize the association between genetic susceptibility to RA and the prespecified disorders, with the results expressed as odds ratios (ORs) and 95% confidence intervals (95% CIs). RESULTS: The wGRS for RA was significantly associated with type 1 diabetes mellitus (DM) (OR 1.10 [95% CI 1.04-1.16]; P = 9.82 × 10-4 ) and multiple sclerosis (OR 0.82 [95% CI 0.77-0.88]; P = 1.73 × 10-8 ), but not with other cardiometabolic phenotypes. In the PheWAS, wGRS was also associated with an increased risk of several autoimmune phenotypes including RA, thyroiditis, and systemic sclerosis, and with a decreased risk of demyelinating disorders. In the IVWR meta-analyses, RA was significantly associated with an increased risk of type 1 DM (P = 1.15 × 10-14 ), with evidence of horizontal pleiotropy (Mendelian Randomization-Egger intercept estimate P = 0.001) likely driven by rs2476601, a PTPN22 variant. The association between type 1 DM and RA remained significant (P = 9.53 × 10-9 ) after excluding rs2476601, with no evidence of horizontal pleiotropy (intercept estimate P = 0.939). RA was also significantly associated with type 2 DM and C-reactive protein levels. These associations were driven by variation in the major histocompatibility complex region. CONCLUSION: This study presents evidence of pleiotropy between the genetic predisposition to RA and associated phenotypes found in other autoimmune and cardiometabolic disorders, including type 1 DM.

Genetic fine mapping of systemic lupus erythematosus MHC associations in Europeans and African Americans.

Genetic variation within the major histocompatibility complex (MHC) contributes substantial risk for systemic lupus erythematosus, but high gene density, extreme polymorphism and extensive linkage disequilibrium (LD) have made fine mapping challenging. To address the problem, we compared two association techniques in two ancestrally diverse populations, African Americans (AAs) and Europeans (EURs). We observed a greater number of Human Leucocyte Antigen (HLA) alleles in AA consistent with the elevated level of recombination in this population. In EUR we observed 50 different A-C-B-DRB1-DQA-DQB multilocus haplotype sequences per hundred individuals; in the AA sample, these multilocus haplotypes were twice as common compared to Europeans. We also observed a strong narrow class II signal in AA as opposed to the long-range LD observed in EUR that includes class I alleles. We performed a Bayesian model choice of the classical HLA alleles and a frequentist analysis that combined both single nucleotide polymorphisms (SNPs) and classical HLA alleles. Both analyses converged on a similar subset of risk HLA alleles: in EUR HLA- B*08:01 + B*18:01 + (DRB1*15:01 frequentist only) + DQA*01:02 + DQB*02:01 + DRB3*02 and in AA HLA-C*17:01 + B*08:01 + DRB1*15:03 + (DQA*01:02 frequentist only) + DQA*02:01 + DQA*05:01+ DQA*05:05 + DQB*03:19 + DQB*02:02. We observed two additional independent SNP associations in both populations: EUR rs146903072 and rs501480; AA rs389883 and rs114118665. The DR2 serotype was best explained by DRB1*15:03 + DQA*01:02 in AA and by DRB1*15:01 + DQA*01:02 in EUR. The DR3 serotype was best explained by DQA*05:01 in AA and by DQB*02:01 in EUR. Despite some differences in underlying HLA allele risk models in EUR and AA, SNP signals across the extended MHC showed remarkable similarity and significant concordance in direction of effect for risk-associated variants.

Whole-exome sequencing uncovers oxidoreductases DHTKD1 and OGDHL as linkers between mitochondrial dysfunction and eosinophilic esophagitis.

Eosinophilic esophagitis (EoE) is an allergic inflammatory esophageal disorder with a complex underlying genetic etiology often associated with other comorbidities. Using whole-exome sequencing (WES) of 63 patients with EoE and 60 unaffected family members and family-based trio analysis, we sought to uncover rare coding variants. WES analysis identified 5 rare, damaging variants in dehydrogenase E1 and transketolase domain-containing 1 (DHTKD1). Rare variant burden analysis revealed an overabundance of putative, potentially damaging DHTKD1 mutations in EoE (P = 0.01). Interestingly, we also identified 7 variants in the DHTKD1 homolog oxoglutarate dehydrogenase-like (OGDHL). Using shRNA-transduced esophageal epithelial cells and/or patient fibroblasts, we further showed that disruption of normal DHTKD1 or OGDHL expression blunts mitochondrial function. Finally, we demonstrated that the loss of DHTKD1 expression increased ROS production and induced the expression of viperin, a gene previously shown to be involved in production of Th2 cytokines in T cells. Viperin had increased expression in esophageal biopsies of EoE patients compared with control individuals and was upregulated by IL-13 in esophageal epithelial cells. These data identify a series of rare genetic variants implicating DHTKD1 and OGDHL in the genetic etiology of EoE and underscore a potential pathogenic role for mitochondrial dysfunction in EoE.

Transcription factors operate across disease loci, with EBNA2 implicated in autoimmunity.

Explaining the genetics of many diseases is challenging because most associations localize to incompletely characterized regulatory regions. Using new computational methods, we show that transcription factors (TFs) occupy multiple loci associated with individual complex genetic disorders. Application to 213 phenotypes and 1,544 TF binding datasets identified 2,264 relationships between hundreds of TFs and 94 phenotypes, including androgen receptor in prostate cancer and GATA3 in breast cancer. Strikingly, nearly half of systemic lupus erythematosus risk loci are occupied by the Epstein-Barr virus EBNA2 protein and many coclustering human TFs, showing gene-environment interaction. Similar EBNA2-anchored associations exist in multiple sclerosis, rheumatoid arthritis, inflammatory bowel disease, type 1 diabetes, juvenile idiopathic arthritis and celiac disease. Instances of allele-dependent DNA binding and downstream effects on gene expression at plausibly causal variants support genetic mechanisms dependent on EBNA2. Our results nominate mechanisms that operate across risk loci within disease phenotypes, suggesting new models for disease origins.

Copb2 is essential for embryogenesis and hypomorphic mutations cause human microcephaly.

Primary microcephaly is a congenital brain malformation characterized by a head circumference less than three standard deviations below the mean for age and sex and results in moderate to severe mental deficiencies and decreased lifespan. We recently studied two children with primary microcephaly in an otherwise unaffected family. Exome sequencing identified an autosomal recessive mutation leading to an amino acid substitution in a WD40 domain of the highly conserved Coatomer Protein Complex, Subunit Beta 2 (COPB2). To study the role of Copb2 in neural development, we utilized genome-editing technology to generate an allelic series in the mouse. Two independent null alleles revealed that Copb2 is essential for early stages of embryogenesis. Mice homozygous for the patient variant (Copb2R254C/R254C) appear to have a grossly normal phenotype, likely due to differences in corticogenesis between the two species. Strikingly, mice heterozygous for the patient mutation and a null allele (Copb2R254C/Zfn) show a severe perinatal phenotype including low neonatal weight, significantly increased apoptosis in the brain, and death within the first week of life. Immunostaining of the Copb2R254C/Zfnbrain revealed a reduction in layer V (CTIP2+) neurons, while the overall cell density of the cortex is unchanged. Moreover, neurospheres derived from animals with Copb2 variants grew less than control. These results identify a general requirement for COPB2 in embryogenesis and a specific role in corticogenesis. We further demonstrate the utility of CRISPR-Cas9 generated mouse models in the study of potential pathogenicity of variants of potential clinical interest.

Genome-wide association study evaluating single-nucleotide polymorphisms and outcomes in patients with advanced stage serous ovarian or primary peritoneal cancer: An NRG Oncology/Gynecologic Oncology Group study.

OBJECTIVE: This study evaluated single nucleotide polymorphisms (SNPs) associated with progression free (PFS) and overall survival (OS) in patients with advanced stage serous EOC. METHODS: Patients enrolled in GOG-172 and 182 who provided specimens for translational research and consent were included. Germline DNA was evaluated with the Illumina's HumanOMNI1-Quad beadchips and scanned using Illumina's iScan optical imaging system. SNPs with allele frequency>0.05 and genotyping rate>0.98 were included. Analysis of SNPs for PFS and OS was done using Cox regression. Statistical significance was determined using Bonferroni corrected p-values with genomic control adjustment. RESULTS: The initial GWAS analysis included 1,124,677 markers in 396 patients. To obtain the final data set, quality control checks were performed and limited to serous tumors and self-identified Caucasian race. In total 636,555 SNPs and 289 patients passed all the filters. The pre-specified statistical level of significance was 7.855e-08. No SNPs met this criteria for PFS or OS, however, two SNPs were close to significance (rs10899426 p-2.144e-08) (rs6256 p-9.774e-07) for PFS and 2 different SNPs were identified (rs295315 p-7.536e-07; rs17693104 p-7.734e-07) which were close to significance for OS. CONCLUSIONS: Using the pre-specified level of significance of 1×10-08, we did not identify any SNPs of statistical significance for OS or PFS, however several were close. The SNP's identified in this GWAS study will require validation and these preliminary findings may lead to identification of novel pathways and biomarkers.

Transancestral mapping and genetic load in systemic lupus erythematosus.

Systemic lupus erythematosus (SLE) is an autoimmune disease with marked gender and ethnic disparities. We report a large transancestral association study of SLE using Immunochip genotype data from 27,574 individuals of European (EA), African (AA) and Hispanic Amerindian (HA) ancestry. We identify 58 distinct non-HLA regions in EA, 9 in AA and 16 in HA (∼50% of these regions have multiple independent associations); these include 24 novel SLE regions (P<5 × 10-8), refined association signals in established regions, extended associations to additional ancestries, and a disentangled complex HLA multigenic effect. The risk allele count (genetic load) exhibits an accelerating pattern of SLE risk, leading us to posit a cumulative hit hypothesis for autoimmune disease. Comparing results across the three ancestries identifies both ancestry-dependent and ancestry-independent contributions to SLE risk. Our results are consistent with the unique and complex histories of the populations sampled, and collectively help clarify the genetic architecture and ethnic disparities in SLE.

Identification of a Sjögren's syndrome susceptibility locus at OAS1 that influences isoform switching, protein expression, and responsiveness to type I interferons.

Sjögren's syndrome (SS) is a common, autoimmune exocrinopathy distinguished by keratoconjunctivitis sicca and xerostomia. Patients frequently develop serious complications including lymphoma, pulmonary dysfunction, neuropathy, vasculitis, and debilitating fatigue. Dysregulation of type I interferon (IFN) pathway is a prominent feature of SS and is correlated with increased autoantibody titers and disease severity. To identify genetic determinants of IFN pathway dysregulation in SS, we performed cis-expression quantitative trait locus (eQTL) analyses focusing on differentially expressed type I IFN-inducible transcripts identified through a transcriptome profiling study. Multiple cis-eQTLs were associated with transcript levels of 2'-5'-oligoadenylate synthetase 1 (OAS1) peaking at rs10774671 (PeQTL = 6.05 × 10-14). Association of rs10774671 with SS susceptibility was identified and confirmed through meta-analysis of two independent cohorts (Pmeta = 2.59 × 10-9; odds ratio = 0.75; 95% confidence interval = 0.66-0.86). The risk allele of rs10774671 shifts splicing of OAS1 from production of the p46 isoform to multiple alternative transcripts, including p42, p48, and p44. We found that the isoforms were differentially expressed within each genotype in controls and patients with and without autoantibodies. Furthermore, our results showed that the three alternatively spliced isoforms lacked translational response to type I IFN stimulation. The p48 and p44 isoforms also had impaired protein expression governed by the 3' end of the transcripts. The SS risk allele of rs10774671 has been shown by others to be associated with reduced OAS1 enzymatic activity and ability to clear viral infections, as well as reduced responsiveness to IFN treatment. Our results establish OAS1 as a risk locus for SS and support a potential role for defective viral clearance due to altered IFN response as a genetic pathophysiological basis of this complex autoimmune disease.

Profound loss of esophageal tissue differentiation in patients with eosinophilic esophagitis.

BACKGROUND: A key question in the allergy field is to understand how tissue-specific disease is manifested. Eosinophilic esophagitis (EoE) is an emerging tissue-specific allergic disease with an unclear pathogenesis. OBJECTIVE: Herein we tested the hypothesis that a defect in tissue-specific esophageal genes is an integral part of EoE pathogenesis. METHODS: We interrogated the pattern of expression of esophagus-specific signature genes derived from the Human Protein Atlas in the EoE transcriptome and in EPC2 esophageal epithelial cells. Western blotting and immunofluorescence were used for evaluating expression of esophageal proteins in biopsy specimens from control subjects and patients with active EoE. Whole-exome sequencing was performed to identify mutations in esophagus-specific genes. RESULTS: We found that approximately 39% of the esophagus-specific transcripts were altered in patients with EoE, with approximately 90% being downregulated. The majority of transcriptional changes observed in esophagus-specific genes were reproduced in vitro in esophageal epithelial cells differentiated in the presence of IL-13. Functional enrichment analysis revealed keratinization and differentiation as the most affected biological processes and identified IL-1 cytokines and serine peptidase inhibitors as the most dysregulated esophagus-specific protein families in patients with EoE. Accordingly, biopsy specimens from patients with EoE evidenced a profound loss of tissue differentiation, decreased expression of keratin 4 (KRT4) and cornulin (CRNN), and increased expression of KRT5 and KRT14. Whole-exome sequencing of 33 unrelated patients with EoE revealed 39 rare mutations in 18 esophagus-specific differentially expressed genes. CONCLUSIONS: A tissue-centered analysis has revealed a profound loss of esophageal tissue differentiation (identity) as an integral and specific part of the pathophysiology of EoE and implicated protease- and IL-1-related activities as putative central pathways in disease pathogenesis.

Genetic architecture distinguishes systemic juvenile idiopathic arthritis from other forms of juvenile idiopathic arthritis: clinical and therapeutic implications.

OBJECTIVES: Juvenile idiopathic arthritis (JIA) is a heterogeneous group of conditions unified by the presence of chronic childhood arthritis without an identifiable cause. Systemic JIA (sJIA) is a rare form of JIA characterised by systemic inflammation. sJIA is distinguished from other forms of JIA by unique clinical features and treatment responses that are similar to autoinflammatory diseases. However, approximately half of children with sJIA develop destructive, long-standing arthritis that appears similar to other forms of JIA. Using genomic approaches, we sought to gain novel insights into the pathophysiology of sJIA and its relationship with other forms of JIA. METHODS: We performed a genome-wide association study of 770 children with sJIA collected in nine countries by the International Childhood Arthritis Genetics Consortium. Single nucleotide polymorphisms were tested for association with sJIA. Weighted genetic risk scores were used to compare the genetic architecture of sJIA with other JIA subtypes. RESULTS: The major histocompatibility complex locus and a locus on chromosome 1 each showed association with sJIA exceeding the threshold for genome-wide significance, while 23 other novel loci were suggestive of association with sJIA. Using a combination of genetic and statistical approaches, we found no evidence of shared genetic architecture between sJIA and other common JIA subtypes. CONCLUSIONS: The lack of shared genetic risk factors between sJIA and other JIA subtypes supports the hypothesis that sJIA is a unique disease process and argues for a different classification framework. Research to improve sJIA therapy should target its unique genetics and specific pathophysiological pathways.

Decreased SMG7 expression associates with lupus-risk variants and elevated antinuclear antibody production.

OBJECTIVES: Following up the systemic lupus erythematosus (SLE) genome-wide association studies (GWAS) identification of NMNAT2 at rs2022013, we fine-mapped its 150 kb flanking regions containing NMNAT2 and SMG7 in a 15 292 case-control multi-ancestry population and tested functions of identified variants. METHODS: We performed genotyping using custom array, imputation by IMPUTE 2.1.2 and allele specific functions using quantitative real-time PCR and luciferase reporter transfections. SLE peripheral blood mononuclear cells (PBMCs) were cultured with small interfering RNAs to measure antinuclear antibody (ANA) and cyto/chemokine levels in supernatants using ELISA. RESULTS: We confirmed association at NMNAT2 in European American (EA) and Amerindian/Hispanic ancestries, and identified independent signal at SMG7 tagged by rs2702178 in EA only (p=2.4×10-8, OR=1.23 (95% CI 1.14 to 1.32)). In complete linkage disequilibrium with rs2702178, rs2275675 in the promoter region robustly associated with SMG7 mRNA levels in multiple expression quantitative trait locus (eQTL) datasets. Its risk allele was dose-dependently associated with decreased SMG7 mRNA levels in PBMCs of 86 patients with SLE and 119 controls (p=1.1×10-3 and 6.8×10-8, respectively) and conferred reduced transcription activity in transfected HEK-293 (human embryonic kidney cell line) and Raji cells (p=0.0035 and 0.0037, respectively). As a critical component in the nonsense-mediated mRNA decay pathway, SMG7 could regulate autoantigens including ribonucleoprotein (RNP) and Smith (Sm). We showed SMG7 mRNA levels in PBMCs correlated inversely with ANA titres of patients with SLE (r=-0.31, p=0.01), and SMG7 knockdown increased levels of ANA IgG and chemokine (C-C motif) ligand 19 in SLE PBMCs (p=2.0×10-5 and 2.0×10-4, respectively). CONCLUSION: We confirmed NMNAT2 and identified independent SMG7 association with SLE. The inverse relationship between levels of the risk allele-associated SMG7 mRNAs and ANA suggested the novel contribution of mRNA surveillance pathway to SLE pathogenesis.